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Abstract 

We revisit XQJ ->• ^ ~t~ (with J = 0, 1,2 and Q = b,c) within the Standard Model (SM). The 
electro-magnetic contributions are given in color-singlet model with non-vanishing lepton masses 
ly-j \ at the leading order of v. Numerically, the branching ratios of XQJ ~^ ^ + ^~ predicted within the 

SM are so small that such decays are barely possible to be detected at future BESIII and SuperB 

experiments. We investigate xw ~^ £ + £~ in Type-II 2HDM with large tan/3, and Xb2 ~ > (- + £~ in 

ft' 

Q_|. the Randall- Sundrum model, to see their chance to be observed in future experiments. 
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I. INTRODUCTION 

The leptonic decay of quarkonium, especially the S-wave triplet J /ip or T, plays an 
very important role in particle physics, either due to its clean experimental signal which 
is commonly used to tag J/ip or T in experiments, or due to its simpleness in theoretical 
calculations which offers an ideal place for precise determinations of the non-perturbative 
Non-relativistic QCD NRQCD) matrix elements |l|. Thus, H( 3 Si) — > £ + £~ has been exten- 



sively studied in past almost four decades |3H7|. 

However, the leptonic decays of C-even quarkonia are of less interests, because they are 
generally suppressed in the SM by both of the electromagnetic loop and huge mass of Z 
joson. The leptonic decay of r] c is investigated by the method of light-cone wave function in 
8| and NRQCD factorization at leading order of typical quark velocity v in [9j. The leptonic 
decay of the C-even and P-wave quarkonia (i.e. 3 Pj quarkonium xqj {J — 0, L2)) has been 
studied with the color singlet model by Kiihn et al about three decades ago [10| . 

In this paper, we will revisit xqj ~~ >* C- + C-~ {J = 0,1,2 and Q = b,c) within the SM. We 



employ the color-singlet model to calculate the decay amplitudes as Kuhn et al did in |10 |. 
We get the decay amplitudes with finite lepton mass m^, which are relevant for the study 
the helicity-suppressed decay xqo ~~ >" £ + £~, and XbJ or higher 3 Pj charmonium excitations 
decays to t + t~. As it should be, our results agree with those obtained in [lOj by setting 
rri£ — > 0. Furthermore, we calculate the electromagnetic loop by utilizing the method of 



regions 111. Il2| . It allows us to relate our results to those calculated in the conventional 



NRQCD factorization for the P-wave quarkonia decays and productions as in [13HL5]. 

Phenomenologically, as expected, the decay widths of xqj —> ^ + ^~ are highly suppressed 
in the SM, so that it is barely possible to measure such decays even in now-days or near- 
future high-luminosity colliders, such as BESIII, LHC and SuperB. However, in an era 
longing for the new physics (NP), the smallness of the branching ratios of such decays in 
the SM can be a virtue. Any experimental discovery of such decays could be an indication 
of the NP. Moreover, the quantum numbers of xqj make the NP effects in xqj ~~ >* £ + £~ 
spin-dependent. We consider two scenarios of extensions of the SM which may enhance 
Xbo — > £ + £~ and Xb2 -^ £ + £~'- one is Type-II two Higgs doublet model (2HDM) with large 



tan/3 16], another is the Randall- Sundrum (RS) model for large extra-dimensions 171. Il8|. 



This paper is organized as follows. In SectJTTl we calculate the decay amplitudes for 



Xqj ~~ >* ^ + ^~ within the SM, compare our results with those obtained in [101 ]. and discuss 
briefly the breakdown and restoration of the NRQCD factorization for xqj ~* ^ + ^~ along 
the way that Beneke and Vernazza did for B — > XcjK in 19]. In Sect. IIH[ we calculate the 
X6o -> £ + t~ in Type-II 2HDM, and Xb2 ->■ Pt~ in the RS model. SectjTV]is devoted for the 
numerical results and some phenomenological discussions. Finally, we summarize our work 
in Sect. EJ 



II. THE DECAY AMPLITUDES FOR X QJ -»■ ^ + ^ WITHIN THE SM 



In the SM, the lowest order Feynman diagrams for xqj — > vl~ are two electro-magnetic 
(EM) box-diagrams and a tree-level Z°-exchange diagram, depicted in FIGJIJ Note that 
only the 3 P\ state can decay into a lepton pair via virtual Z° at tree-level, and the tree-level 
neutral higgs exchange diagram for xqo —> ^ + ^~ is neglected. 
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FIG. 1: The Feynman diagrams for s Pj state QQ direct annihi 
the Feynman diagrams in this paper are created by JaxoDraw 



ation to lepton pair in the SM. All 

20]. 



We perform the calculations in the rest frame of quarkonium. The the four-momenta of 
the quark, anti-quark, lepton and anti-lepton are denoted as 



Pi 
Qi 



(£,q), p 2 ■■ 
(^,qi), 92 






(la) 
(lb) 



where E = M XQJ /2 with M XQJ being mass of Xqj- We set 2P as the total momentum of 
QQ, and 2q the relative momentum of QQ, i.e. 



r, P1+P2 , „ n x P1-P2 (n s 

p = — g — = ( E '°)' Q = — 2 — = (°' q )' 
The relative velocity between the heavy quark and anti-quark is defined as v = |q|/mg. 



(2) 



As in 



2ll |22| . we describe briefly the procedure to calculate the P- wave quarkonium 



involved processes within the color singlet model developed by Kiihn et al in [10 1 



To project the spin-triplet part of the amplitude QQ( 3 Pj) — > £ + £ , we use the covariant 
projection operator of spin-triplets [23[ 

ir(spin-triplet) = - -— 1— - - (f»+^ + m Q )(# + £)t"(?- tf - ™q) , (3) 

to replace the quark-anti-quark spinor bilinear, where tuq is the heavy quark mass. We 
expand the spin-triplet amplitude as series of the small momentum q 

iM u (q) = iM»(0) + q^iM u (q)\ q=0 + ■ ■ ■ • (4) 

Then, the amplitude at the leading order of v for a 3 Pj state decay to £ + £~ can be obtained 
through 

iMfPj -4 £ + r} = \q\U^( 3 Pj)-^iM v (q)\ q=0 , (5) 



with the projectors 



n„( 3 Po) = ^(-^ + ^), (6a) 



^3r % 



n ^ ( l} = V2E e ^ PPe ^ (6b) 

n^( 3 P 2 ) = e^, (6c) 

where e CT and e M " are respectively the polarization vector and tensor for the 3 P± and 3 P2 
states, which satisfy the relations P M e^ = 0, P^e^ u = and t^ = 0. The polarization 
summations are 10J 

X)c^ = -^ + ^ = P^, (7a) 

/ J CfivCaP = ~ [P^aP^/3 + ^fi^ua] ~ Ai/Pa^ • (7b) 

After convoluting the amplitudes in flHJ) with the radial wave-function for %qj in momentum- 
space, we reach the final amplitudes for xqj ~ > ^ + ^~ which is 

d_ 



tM[ XQ j^£ + n = -\hM XQJ N c \l—RXm^Pj)—iM»{q)\ q=0l (8) 



where the pre-factor J2M XQ1 originates from the relativistic normalization of the xqj state, 
^/W c from summation of the colors with N c = 3 being the number of colors, and -R'(O) is 



derivative of the radial wave function of xqj at origin in coordinate-space [10]. 



A. Calculation of the box diagrams 



We ca' 



to 
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culate the two box- diagrams depicted in FIG. [T]by use of the method of regions 
12| . The leading regions in these box diagrams for the 3 Pj state decays into lepton 
pair in expansion of the relative velocity v, are: 1) hard region where each component of 
the momenta of both photons in the loop at order of tuq] 2) ultra-soft region where each 
component of the momentum of one photon at order of ttiqv 2 . In each region, the power 
counting rules of all momenta are clear so that we can perform the f-expansion for the 
loop-integrands straightforwardly Then we integrate the expanded loop-integrand over the 
whole momentum space to get the contributions from each region, which sum reproduce the 
complete amplitudes in a series of v . 

With all the techniques described above, after some tedious but straightforward calcula- 
tion, we get the amplitudes from the hard region which are infrared (IR) divergent. Here we 
use the dimensional regularization (DR) to regulate the IR divergence. The corresponding 



amplitudes are 



iM hard [ 3 P 



2 v / 6a 2 e| 
3(3m 3 Q 



In 



1 + 13 2/3 - In i±f x i 



P 



ft 2 



x \q\m e u(qi)v(q 2 ) 



h finite terms 

cm 



(9a) 



iM hatd \ 3 P 1 ^e + i 



2^4 

/3m 2 Q 



In 



fi , 2/?-lni±f 



1-/3 



+ 



il 2 



x IquteihsMfc 



h finite terms 

em 



(9b) 



iM^ d [ 3 P 2 ^£ + t 



2V2a 2 e 2 



(3m 3 



(> 



In 



1+0 2/3-lnf±|^ | 



P 



+ 



J 2 



x h\u(qih a v(q 2 )q ^e 



em 



finite terms 



(9c) 



where (3 = J 1 — m 2 /m 2 Q is the velocity of lepton, a = e 2 /(4ii) the fine-structure constant, 
and eg the electric charge in unit of elementary charge e for the heavy quark Q. We do not list 
the explicit expressions for the finite terms above, since they are somewhat scheme-dependent 
and therefore meaningless individually. One should notice that «A^hard[ 3 -Po - > £ + £~] is 
proportional to rri£ from the helicity suppression. 



In a similar way, we get the contributions from the ultra-soft region to the amplitudes 
which are ultraviolet (UV) divergent. Here we also use the (DR) to regulate the UV diver- 
gence. The corresponding amplitudes are 



iM us [ 3 P -)• 



2 v / 6a 2 e| 
3(3m 3 Q 



In 



1 + 2/3-lni±f\ i 



1-/3 



+ 



/? 2 



euv 



+ finte terms 



iMj?P x -)• 



x |q|m £ M(gi)f(g 2 ) 
2a 2 e 2 r 



(10a) 



/3?7l/ 



In H — — 1 h finte terms 



1-/3 



/^ 2 



euv 



X |q|M(?l)75^(?2) 



(10b) 



iMJ 3 P 2 -» 



; 2yg^4 

/3m 3 



In 



l + /3 2/3-ln±±f\ i 



1-/3 



+ 



/3 2 



euv 



finte terms 



x Hu(qih a v(q 2 )q p e 



a/3 



(10c) 



It is easy to see from (j§j) and (TTUl) . that the divergent parts of hard and ultra-soft parts have 
opposite signs. Thus, the whole amplitudes for 3 Pj — ¥ £ + £~ due to the EM interactions 



iM em [ 3 Pj -»■ 



^^^hard[ 3 i ^ J "► ^ + r] + zA^ us [ 3 Pj 



(11) 



are finite. 



In all, we have 



iMemlXQJ ~* 




e^a 2 



-i,/2M XQJ N cX l-R\0)^fjx{ 



meu(qi)v(q 2 ) , J = 0, 

mQu(qi)^v(q 2 )^ k v(q 2 ) , J = 1 , 
%i)7^fe)g 1/3 e^ , J = 2. 

(12) 

The finite coefficients /j( J = 0, 1, 2) are from the sums of ([9]) and ( TTUl) . their explicit analytic 
expressions are 

2v/6 (2/3 - (1 - /3 2 ) In (f±|)) In (^) v^ ((/3 2 + 2) Li 2 (g) - 2 (/3 2 - 1) Li 2 (^f 



/, 



+ 



+ 



2 v / 6/3 ; 



3/3 3 



1-/3; 6^/3^ 



3/3 3 
[tt 2 (5/3 2 - 2) + 12/3 ln(l - /3 2 ) + 24/3] 



7p ln (1=7) [( ~ 2 ~ 3ot)/32 ~ 4(1 " ^ ln(4/3) + 2(1 " ^ 2) ln(1 " f) + 2] 



(13a) 



./l 



2 (2/3-(l-/3 2 )ln(fii) 



2(l-/3 2 
/? 3 



/3 3 

In 



H7) 



(I-/? 2 ) 



/3 3 



Li 2 



1-/3 

/3 + 1 



4/3 

/3 + 1 



- 1 In 



/3 + 1 
1-/3 



7T _ 

12 



(13b) 



^^f^-d.^/^ 



K^) 



+ 



/5 3 V Vi-/3, 

V2 (/3 2 - I) 2 Li 2 (gg) 2y^ (/3 2 - 1) Li 2 (jg ) (1 - 3/3 2 ) (1 - /3 2 ) In 2 (^f) 

^ 272> 

+ (/3 4 - l)wr + 4/3 4 



1 



In 



/3 5 
1-/3 



2/3 2 (l-/3 2 ) 41n2-l-ln 



l-/3 ; 
/3 2 



\/2>~V/3 + l 
+ ^r [-tt 2 (1 - /3 2 )(1 + 3/3 2 ) + 4/3(3 + /3 2 )(2 In 2 - «r) 

+16/3 3 - 12/3(1 + /3 2 )ln(l-/3 2 



(13c) 

Here, uj = M XQJ — 2m<g + ie is the binding energy of Xqj> an d Li 2 (x) = — /j" hi(l — x)/x is 
the Spencer function . 

In the massless lepton limit r = mj/mQ = 1 — /3 2 — > 0, 



/o-> 



161n(— ) +31n 2 r + (4 + 6«7r- 121n2)lnr 



/i->41n 



2^6 

+ 7T 2 + 8 + 12 In 2 2 - 12m In 2] 



Li,' 



/2 _,_^f ln2 _31n^-.l + i ff 

3 V u; 



(14a) 
(14b) 

(14c) 



One can see that /o diverges when r — > while /i )2 remain finite in the same limit. However, 
the decay amplitude for xqo ~ > ^ + ^~ still vanishes when r — > 0, since we have singled out 
the helicity-suppression factor mi. Finally, we reproduce the results of Kiihn et al in [lOj: 



iMemiXQo -► ^ 



->-o, 



15a) 




2 ^,2 



4e^a 



™>q. 



> - h /2M XQ1 /V 4 /- J R'(0)^^ln-^w(g 1 )75^(g 2 ), (15b) 

(15c) 
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FIG. 2: The real and imaginary parts of /j as functions of r = rrij/rri^. 

Taking Xcj as examples, we plot the real and imaginary parts of the coefficients fj 
(J = 0, 1, 2) as functions of r = mj/m^ in FIG. El We roughly take the m c = 1.65 GeV and 
M Xcl23 = 3.41,3.51,3.56 GeV. /i^ are almost flat when r < 10~ 2 , but decrease evidently 
where r is greater than several percent. This means that our results is useful to make more 
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accurate predictions for Xbj or higher radial excitations of Xcj decays to t + t~ . 

Moreover, theoretically, the imaginary parts of fj originate from the on-shell immediate 
states. In FIG. HJ we see that the imaginary part of f\ vanishes because a massive (axial)- 



vector cannot decay to two photons due to Yang's theorem 24( . However, one should notice 
that, if the binding energy u can be negative, the lnrng/w term in f\ can contribute an 
imaginary part. And the In a; term comes only from the ultra-soft part in ( fit)]) . As Kuhn et 



al have pointed out in [10|], that this effect is related to the El transition Pj — > Si + 7. 

B. The neutral current contributions 

For completeness, we present the neutral current weak interaction contributions to \Qi ~~ >* 
£ + £~ as follows. 



iM weak [ X Qi -»• f + r] = ± l G F J^-R'(0)u(q 1 ){ (1 - 4sin 2 ^ - 75) v(q 2 ) , (16) 

V nM ^ 

where the pre-factor "+" is for Q = c, "— " for Q = b, Gp the Fermi constant, and 9w 
the Weinberg angle. As we will see in the numerical analysis, the neutral-current contribu- 
tion may play an important role, especially in Xbj —*■ £ + £~ since its electromagnetic decay 
amplitude is further suppressed by e 2 = 1/9. 

In all, the decay widths for xqj ~^ £ + £~ within the SM can be written as 



t[xqo -► £ + n = A \ Q , /? 3 (i - /3 2 )i/ | 2 |it;'(o)i 2 , (17) 



47r 2 m 



nx Q i-+t + n = ^-i\R'm 



4 4 

a 

4 
> 



R lD//m|2 
4V2e 2 Q a 2 G F 



7T 2 



, |/ 1 | 2 /3 2 + G|(2« + (3-/3 2 )^) 

m Q 



rri^ 



P 2 9aMh 



(18) 



N n ef,a 4 .. ., , 



r[XQ2^£ + n = —^^(5 -2/3 2 )|/ 2 | 2 |i?'(0)| 2 , (19) 

ZXMX Tfhr\ 

where g v — 1 — 4 sin 2 9w, g a = — 1, arid "±" correspond to Q = c and Q = b, respectively. 

C. Connections to the NRQCD factorization 

Before we get into phenomenological applications of our results obtained above, we would 
like to translate our calculations into the language of the effective field theory, and see what 



we can get from such comparisons. 

Heavy quarkonium decays involve three well-separated intrinsic scales tjiq, rriQV, ttiqv 2 . 
The NRQCD is a suitable and powerful effective field theory to describe the heavy quarkonia 
production and decays [l|. By integrating out the hard fluctuations around the scale m,Q, 
the effective Lagrangian for the leptonic decay of a non-relativistically moving heavy quark 
pair QQ can be written as 

5C =^0( 3 5 1 ) + J2 l ^T lo ( ip j) + ■ ■ ■ > (2°) 

y J=l,2,3 v 

where f( 2S+1 Lj) are the short- distance coefficients supposed to be finite, and the effective 
operators are 

(21a) 
(21b) 

(21c) 

(21d) 

Here, following the conventions adopted in [19], we employ the four-component spinors fields 
i[) and x to represent the the non- relativist ic heavy quark and anti-quark, respectively. These 
fields satisfy 7°^ = ip and 7°x = — Xi an d are equivalent to the conventional NRQCD two- 
component fields used in |l|. The Latin indices i,j run over the spacial indices 1,2,3, and 
(i,j) means the traceless part of a symmetric tensor. And 

X^DT^^x^DT^-x^DTil;, (22) 

IdTt^edTe-idTt, (23) 



0( 3 ^i) = xh 


tyhA 






o(?Po)~ 2 ™V( \V)*U 


0{*P 1 ) = -±= X ] (~ D l [7<,7ib B ) W 




o( 3 p 2 ) = x f 


(4)<3V 





where V denotes any Dirac structure, D l = d l — ig s A a,l T a the covariant derivative in QCD. 
Note that the P-wave operators 0( 3 Pj) are v -suppressed relative to the S"-wave operator 
0( 3 Si) due to the power-counting rules in NRQCD. 

Neglecting the contributions from weak interaction, at the lowest order of the strong 
coupling a s , we have the expectation /( 3 5i) at 0(a) while f( 3 Pj) at 0(a 2 ). Naively, one 
would expect that 

^M[x Q .J -> ttT] = i^^-{£ + r\0( 3 Pj)\ XQ j) , (24) 

m Q 
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within the frame of NRQCD at the leading order of v. Then, the validation of (l24j) implies 
a factorization, in which the short- distance contributions are absorbed into f( 3 Pj) while all 
the long-distance contributions are absorbed into the matrix-elements of 0( 3 Pj). 

However, one should notice that, within the NRQCD, an ultra-soft photon can interact 
with the quarks and leptons, and such interaction can be described by 25| 



5C US = ee Q ^(x) [A° cm (t) + x • E(t)] i/>(x) 

-ee QX \x) [A° cm (t) + x • E(t)] X (x) , (25) 

where A® m is the electro-potential, and E the electro-field strength. Both of A° cra and E 
have been multipole expanded. In real calculations of the amplitudes, the A term in (1231) 
can be dropped out either by choosing the temporal gauge A® m = or by the automatic 
cancellations in the amplitudes in other gauge choice. The E term in (1251) is actually the 
electro-dipole interaction, which is at the order 0(v) and conserves the spin but change the 
orbit angular momentum L by one unit. It implies that the ultra-soft photon interaction in 
(|2"5"|) can transform a 3 Si state into a 3 Pj state with a price of 0{y) suppression. 

Therefore, the correct decay amplitude for xqj ~~ >* ^ + ^~ at the lowest order of v within 
the NRQCD should be written as 



J ^(t + nO( 3 Pj)\XQj)U + i f -^ ,~ . -v— v./i^oop 



(20) 



The matrix-element of 0( 3 S\) above can be depicted by the NRQCD Feynman diagrams in 
FIG. 

One familiar with the method of regions, can immediately recognize the matching equa- 
tions 



J( 3 P 



■LLi f)+(}-\fi(3r>\\,, / „ \\ — ; \a. .k, v e+ 



mil 



\0( 3 Pj)\ XQ j)\tree = tM havd [ X QJ ~> t T] , (27a) 



./W (£n-\0( 3 S 1 )\ XQ j)\^ loop = iM us [Xqj -»• Ft-]. (27b) 

m Q 

Consequently, one can find that the "short-distance" f( 3 Pj) does contain IR divergences 
which breaks down the naive factorization in ( 12^1) . 

Of course, the breakdown of the conventional NRQCD factorization for many P-wave 
quarkonium involved processes is not new. Taking B — > XcjK for instance, the QCD 
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Q 

Q 





r 



FIG. 3: The ultra-soft loop corrections to (0( 3 Si)}. 



factorization breaks down 26H31|. However, Beneke and Vernazza showed in 19], that the 
factorization for B — > XcjK can be restored, by considering the contribution from S"-wave 
color-octet operators in which the chromo-El transition analogue to ([2"o]) plays a crucial role. 
Therefore, the identifications in ( 1271) and fmiteness of (1261) can be regarded as an appli- 
cation of the idea developed by Beneke and Vernazza, to restore the factorization. However, 
since it deviates from our final phenomenological goal, we would like to stop the further 
discussions along this line. 



III. POSSIBLE IMPACTS FROM NEW PHYSICS BEYOND THE SM 



As we have seen, xqj ~ > ^ + ^~ is highly suppressed in the SM, due to either the EM loop 
or the large mass of Z°. For xqo —> £ + £~, it suffers more suppressions in the SM due to 
the helicity selection rules. The tininess of the branching ratios make such decays sensitive 
to the possible new physics beyond the SM. If the quantum numbers of the new particles 
in the SM extensions match those of XQJi an d the couplings among them are enhanced in 
some way, we may have a chance to find the hints of new physics in XQJ -.«- In this 
section, we consider two kinds of models: Type-II 2HDM with large tan p [16] , and the RS 
model of the large extra-dimension 13, [18]. 
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A. xqo -»■ £ + t~ in Type-II 2HDM 

Type-II 2HDM is one of the most studied extensions of the SM. It shares almost the 
same higgs sector interactions with the minimal super-symmetric SM (MSSM). The general 
Yukawa couplings among the fermions and the lightest neutral higgs h in 2HDM can be 
written as 

-2HDM 9 i ( COSa rTTT sma 7=.7-. Sma 7I, 



Cy = ~7^ h rn v —-UU - m D -DD - m e -U , (28) 

2m\y \ sinp cosp cosp / 

where U denotes for the up-type quarks, D for the down-type quarks, £ for the charged 
leptons, g for the SU(2) L gauge coupling, a for the mixing-angle of the neutral higgs, and 
tan (3 = Vd/v u with Vd, u being the vacuum expectation values of two higgs doublets coupled 
to the down-type and up-type quarks respectively. 

Straightforwardly, we have the amplitudes for xqo ~~ >* £ + £ + via h in Type-II 2HDM are 



6N C sm2a m c mf,_ I 3 

7^ Gf sin 2pm 2 U{qi)V ^ I "V^- 



tM 2HDM [Xco -»- tr] = i^\^G /ZlS 2 e u(qi)v(q 2 ) ( -\lizlt{0) ) , (29) 



QN C sin 2 a m b m e _ ( /3 

^ Gf cos* (3ml u{qiHq2) Hto- 



tM 2HDM [ Xb o -► £ + H = -i^l^G F _ 9 a * *&)<<*) I -\lTZ#(P) I > (30) 



where m^ is the mass of the lightest neutral higgs. Here we have used the relation 4a/2GV = 

9 2 /m 2 w 

In the large tan/3 scenario of Type II 2HDM, i.e. tan (3 » 1 and a ~ /?, the amplitude for 

Xw ~~ >* £ + £~ is enhanced by the factor m^m^ tan 2 j3, which may compensate the suppression 

from the factor m\, while Xco —^ £ + £~ does not receive such enhancement. Thus, in the 

numerical analysis below, we will consider only xw -* £ + £~ in Type-II 2HDM with large 

tan/3. 



B. XQ2 ->• t + £~ in the RS model 

As a potential solution to the hierarchy problem, the Randall-Sundrum(RS) model [17| 
predict that a TeV Kluza-Klein (KK) resonances may couple to the SM particles. The 



corresponding effective Lagrangian is 18] 



H 



-int — V ^SMi (31) 
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where h^ is the KK graviton field, Tsm the SM energy-momentum stress tensor, 



K 



k/Mpi the effective coupling constant with Mpi = Mpij \/&k the reduced Plank scale, 
k ~ Mpi is the space-time curvature in the extra dimension, mo = ke~ klTrc is a mass scale at 
the order of TeV, and the KK graviton masses are rrii = moXi, where r c is the compatification 
radius of the extra dimension, and Xi are roots of Bessel function J\(x). 

The stress tensor for fermions is proportional to ip (— |) D ^7^ ip, which matches to the 
quantum number of a 2 ++ quarkonium. Thus, XQ2 ~^ d- + d-~ can happen at tree-level. Here 
we consider the contributions from the lowest KK excitations of graviton. A straightforward 
calculation shows that 

iM RS [XQ2 -> £ + H =z2V2^u( qi )^v(q 2 ) {qi ~ 2 q2h e a , x (-JE#(p)\ , (32) 



and the resulted decay width is 

l#(0)r(-V^(5-2/3 2 ), (33) 

\mij m,Q 



^ s _32/3 3 iV C| .^ |2 / K \ 4 m 4 Q 



5tt 2 

where mi is the mass of the lightest RS graviton. Though the RS graviton contribution is 
TeV scale suppressed, nevertheless, XQ2 — > C- + C-~ may be sizable if k is not too small. Since 
(133]) indicates that T RS is proportional to rrin, we only consider Xb2 -^ ^ + ^~ in numerical 
analysis below. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

Here we present the numerical results based on our calculations of xqj ~^ £ + £~ within 
the SM and beyond. We take the following values for the fine structure constant, Fermi 
constant and Weinberg angle 

a = 1/137, G F = 1.16 x 1(T 5 GeV, sin 2 B w = 0.231, (34) 

in our numerical analysis. 
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A. XQJ -> £ + £~ within the SM 



1- Xcj 



For XcJ, w e take 



M(x c o) = 3.41 GeV, M{ Xc i) = 3.51 GeV, M( Xc2 ) = 3.56 GeV. 
r(x c o) = 10.3 MeV , r(xci) = 0.86 MeV , L( Xc2 ) = 1.97 MeV . 



(35) 
(36) 



The binding energies are taken as u = M(x c j) — 2m c . The derivative of Xcj radial wave 
function at the origin is taken as |i^ c (0)| 2 = 0.050GeV 5 [191 ]. 

In TABLE [U we show the branching ratios of Xcj — > e + e~ and /i + /i~ with m c = 1.65GeV. 
One can see little difference between the Br[x c i,2 — > e + e~] and Br[x c i,2 — * ^ + ^}i but large 
differences between Br[x c0 — y e + e~] and Br[x c0 — > n + n~], because Xco —> £ + ^~ is helicity 
supressed. 

TABLE I: The branching ratios of Xcj -^ i + t~ . 



Decay Channels 


QED 


QED+weak 


Xco -> e + e" 


3.19 x 1(T 13 




Xco ->• V> + fi>~ 


7.06 x KT 10 




Xc\ -> e + e" 


4.55 x 10~ 8 


3.67 x 1(T 8 


Xc\ -> M + M~ 


4.44 x KT 8 


3.60 x 10~ 8 


Xc2 -> e+e" 


1.37 x 10~ 8 




Xc2 ->• /^ + /i~ 


1.32 x KT 8 





To show the dependence of our results on the charm quark mass, we also plot Br[x c j — y 
e + e~] as functions of m c in FIGJU The curves show the strong dependence of the branching 
ratios on m c , or more precisely the binding energy u, which is mainly due to the lnm c /w term 
in the decay amplitudes. And the cusps of the curves also originate from the discontinuity 
of this logarithmic term. 

The most of Xcj events produced at BEPC-II are from the radiative decays ip(2S) —y Xcj + 
7. BEPC-II will produce about 3.2 x 10 9 ip{2S) per year. With Br[ip(2S) -)- Xco+l] = 9.62%, 
Bt[iP(2S) -)■ Xd + l} = 9.2% and Br[^(2S) -> X c2 + 1] = 8.74%, we expect that there will 
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FIG. 4: Br[x 



cJ 



e ' e 



as functions of m c . 



be about 2.89 x 10 8 Xco, 2.76 x 10 8 Xei and 2.62 x 10 8 \c2 produced at BEPC-II annually. 
From the results listed in TABLE 0] it seems that we have no possibility to see Xco —*■ £ + P~, 
but marginal possibility to see x&,2 ~ > ^ + ^~ at the BES-III experiments. 

2. Xbj^t + l- 



For the Xbj leptonic decay, we take 



M X60 = 9.86 GeV, M Xbl 



9.89 GeV, M. 



:vw> 



9.91 GeV. 



(37) 



and |i?'(0) X67 | 2 = 1.42 GeV 5 given by the Buchmuller-Tye potential model [21|. The total 
decay widths of the Xbj are absent in PDG review, we list the decay widths of Xbj —* ^ + ^~ 
in TABLE |TT] by setting m b = 4.67 GeV. 

One can see that the decay widths of Xm.,2 -* t + t~ are significantly different from those of 
Xbi,2 —> e + e~ and /i + /i~. This makes our recalculations of the EM box diagrams meaningful. 
One can also find that the Z°-exchange dominates Xbi ~~ >* C- + C-~ ■ 

To see the uncertainties of our results due to the bottom quark mass m&, we also list the 
decay widths for Xbj —$■ ^ + ^~ in TABLE [TTT] with several different values of m b . 



16 



TABLE II: The decay width of Xbj -»• i + t~ in units of GeV. 



Decay Channels 



QED 



QED+weak 



Xbo 
Xbo 
Xbo 



e T e 



/*■> 



T + T 



1.69 x 

4.98 x 
1.11 x 



10 



10 



-17 



-14 



10 



-12 



Xbi 
Xbi 
Xbi 



e ' e 



V V 



T + T 



1.16 x 
1.16 x 
7.20 x 



10 



-12 



10 



-12 



10 



-13 



2.19 x 10" 11 
2.18 x 10 -11 
1.66 x IO -11 



Xb2 

Xb2 

Xb-2 



e ' e 



fi^fi 



T + T 



9.08 x 
9.02 x 
5.52 x 



io- 13 
io- 13 
io- 13 



TABLE III: The decay width of Xbj ~ 


■> l + d. in units of GeV with different values of m&. 


m 6 /GeV 4.6 


4.8 5.0 



Xbo 
Xbo 
Xbo 



e ' e 



W 



T + T 



1.83 x IO" 17 
5.28 x 10" 14 
1.12 x IO" 12 
2.09 x 10" 11 
2.09 x 10" 11 
1.58 x IO" 11 
8.06 x IO" 13 
8.00 x IO" 13 
4.80 x 10" 13 



1.47 x IO" 17 
4.57 x IO" 14 
1.13 x IO" 12 



1.22 x IO" 17 
4.54 x 10~ 14 
1.70 x IO" 12 
3.03 x IO" 11 
3.03 x IO" 11 
2.33 x IO -11 
4.00 x IO" 12 
3.98 x IO" 12 
2.68 x 10~ 12 



Xbi 
Xbi 
Xbi 



e ' e 



H^fi 



T + T 



2.45 x IO" 11 
2.45 x IO" 11 
1.87 x IO" 11 



Xb2 
Xb2 
Xb2 



fi + fT 



T + T 



1.25 x IO" 12 
1.24 x IO" 12 
7.89 x IO" 13 



B. Impacts from new physics 

We now consider the numerical impacts on xqj ~ > £ + £~ from the new physics effects. 
A complete analysis with considering all possible parameters and constrains to parameter 
space seems rather complicated and deviate from the main topic of the paper. We will 
perform our analysis by choosing some specific values for the new physics parameters. 
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1- Xbo^ t + t~~ in Type-II 2HDM 



In general type-II 2HDM, we have to consider three individual parameters: rrih, tan/3 
and a. Recently, the higgs mass has been excluded in a broad range of mass parameters. 
Recently, both the CMS and ATLAS collaborations observed an excess of higgs-like events 
with 2-3 sigma around 125 GeV in pp collisions at the Large Hadron Collider (LHC) at 



t/s = 7 TeV |32|, |33|. In the following analysis, we roughly take rrih around 125 GeV, 



although the experimentally observed events may not be referred to the lightest neutral 
higgs in 2HDM. 



Decay width/GeV 
6.x 10 




FIG. 5: The decay width as functions of r = m 2 /m 2 with tan/3 = 10 and different choices of rrih- 



We first take a ~ /3, while the mass of neutral higgs rrih arid tan /3 are free parameters. 
We plot the decay width of xw — > £ + £~ as a function of r = mf/mf with tan/3 = 10 and 
different values of rrih, in FIGJ5J and also plot the decay width as a function of rrih with 

tan/3 = 5,10,15 in FIG. H 

It is easy to see that, when tan /3 = 10, the decay width is much larger than the SM results. 
If rrih is taken to be 125GeV, the decay width of Xbo —> t + t~ will be about 5.07 x 10 _12 GeV, 
which is about three times of the SM prediction 1.58 x 10~ 12 GeV. Thus, if Type-II 2HDM 
is a true theory of our world, we may have better chance to observe Xbo — > £ + £~ ■ 

In general 2HDM, the neutral higgs mixing angle a is not strongly correlated with /3. To 
see how the decay width of xw -^ £ + £~ depends on a, we plot Xbo -^ t + t~ decay width as 
a function of a with tan/3 = 5, 10, 15 respectively at m h = 125GeV in FIGJ3 When a = 0, 
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FIG. 6: Decay width as functions of m^ with different choices of tan/3. 

the decay width reaches its minimum value 1.58 x lCT 12 GeV, which is just the SM result, 
and it reaches its maximum value when a — (3, which is just the case we have discussed 
above. 



Decay w 


idth/GeV 


tan/3=5 


7.0 xur 12 






5.0 xur 12 




tan6=10 




\ 


tan8=15 


3.0 xur 12 


\ 






FIG. 7: Decay width as functions of a with different choices of tan/3. 



-2- Xb2 -^ f + £ in the RS model 

In the RS model, the most important parameters are k = k/Mpi. Recently, the AT- 
LAS collaboration has reported the 95% CL. lower limit on the mass of RS graviton for 
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various values of k/M Pl , which are 0.71, 1.03,1. 33, 1.63TeV for k/M Pl = 0.01,0.03,0.05,0.1, 
respectively J34|. 

Here, we only consider the decay width induced by the lightest KK excitation of the 
graviton with neglecting the interference between graviton-exchange and SM contribution. It 
will not prevent us getting a qualitative observation. We plot the decay width of Xt2 —* £ + ^~ 
via lightest KK graviton as a function of m\, with k/Mpi = 0.01,0.03,0.05,0.1 in FIG]8j 
while the lepton mass me is taken to be 0. From FIG. [HJ TABLE [TT] and |III|, one can see 
that, even for k/Mpi ~ 0.1 and vti\ ~ 200 GeV, the KK graviton exchange cannot compete 
with the SM contributions in Xb2 —*■ i + £~- Thus, we have no chance to look for any hints of 
the RS model in such decay. 



Decay width/GeV 

k 

io- 1 




k/Mpi=0.01 



k/Mpi=0.03 



k/M P ,=0.05 



500 



1000 



1500 



2000 



,m,/GeV 



FIG. 8: The decay width of Xb2 ~ > ^ + ^ as function of m\ with different choices of k/Mp\. 



V. SUMMARY 

In this paper, we recalculate xqj ~~ >* £ + £~ within the SM by considering the finite mass of 
the leptons. The suppression of such decays in the SM make them sensitive to the NP. In the 
future experiments where massive quarkonia are produced, such rare decays of quarkonia 
could be a new play-ground of NP hunters other than high-energy colliders or super flavor 
factories. We investigate Xbo —* ^ + ^~ m Type-II 2HDM, and Xb2 — >• ^ + ^~ in the RS model. 
We find that in the large tan/3 limit, we may have better chance to observe Xbo -^ t + t~ 
in Type-II 2HDM than in the SM, and no chance to find the hints of the RS model in 
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Xb2 —> t^^ ■ It could be interesting to investigate, in which kind of extensions of the SM, 
the leptonic decays of xqj might be enhanced so that we can observe them. 
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